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1 Introduction 

Starting in sixties an active group of physicists under the guidance of Prof. 
K.A.Ter-Martirosyan was creating the theory of high-energy processes in 
QCD. From the beginning the key element of this theory is the notion of 
Regge trajectories and in particular of the pomeron trajectory, which have 
been introduced phenomenologically. In this talk I review the problem of 
spectrum and Regge trajectories as it can be derived from the nonperturba- 
tive QCD dynamics. 

The problem of the QCD spectrum is the central issue in the nonper- 
turbative QCD and is intimately connected to problems of confinement and 
mass generation in QCD. These latter issue make QCD so much different 
from QED, in addition explicit form of Nonperturbative Interaction (NPI) 
was not known for a long time. Recently with the introduction of the Field 
Correlator Method (FCM) [Q, |^ the situation has changed favorably, since 
the NPI is defined there in a gauge-invariant way and the simple field cor- 
relators, which are sufficient for all dynamical calculations are known from 
lattice data [§, ^ and analytic results p. 

In view of all this one can put the problem of the QCD spectrum in the 
most general framework, as the many-channel problem of bound states of 
quarks, antiquarks and valence gluons with all possible mixings, the states 
being stable in the limit oo, and acquiring decay widths when Nc = 3. 
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It is the purpose of this talk to describe the QCD spectrum in the Hamilto- 
nian language, using FCM and the Nc = oo limit. In addition FCM provides 
another convenient limit - the gluon correlation length Tg tending to zero, 
while the string tension a is kept constant. As a result the universal Hamil- 
tonian becomes local and has a transparent structure for any number of 
constituents. As will be seen, it contains only two fixed input parametersicx 
and as (or Aqcd) in addition to current quark masses and is able to predict 
any meson, baryon, hybrid and glueball states and their mixings (for a review 
see 0). 

In doing so one solves the following puzzles: 1) constituent mass of quarks 
and gluons 2) meson Regge trajectories with a correct slope and intercept 3) 
the explicit notion of the valence gluon 4) hybrid spectra and Regge trajec- 
tories. 

Remarkably all calculated spectra are in very good agreement with lat- 
tice data and experiment which gives an additional justification for FCM. 
However immediate consequences of the whole approach are far more reach- 
ing. They include a new formulation of the perturbation theory, the so- 
called Background Perturbation theory (BPTh) with saturating at 
small Euclidean momenta instead of diverging near the Landau ghost pole. 
The whole structure of QCD becomes interconnected with the spectrum via 
the quark-hadron duality and hybrids and glueballs play a very important 
role in scattering and decay. 



2 Hamiltonian 

There are two possible approaches to incorporating nonperturbative field 
correlators in the quark-antiquark (or 3q) dynamics. The first has to deal 
with the effective nonlocal quark Lagrangian containing field correlators 
From this one obtains first-order Dirac-type integro-differential equations 



for heavy-light mesons p], light mesons ||T0| and baryons |T^, |TT|. these 
equations contain the effect of chiral symmetry breaking which is directly 
connected to confinement. 

The second approach is based on the effective Hamiltonian for any gauge- 
invariant quark-gluon system. In the limit Tg ^ this Hamiltoniam is simple 
and local, and in most cases when spin interaction can be considered as a 
perturbation one obtains results for the spectra in an analytic form, which 
is transparent. 
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For this reason we choose below the second, Hamiltonian approach \T}, 
r3[ . We start with the exact Fock-Feynman-Schwinger Representation for 
the qq Green's function (for a review see [0), taking for simphcity nonzero 
flavor case 

G\7 = I dsij^ ds2izU{Dz),yr-^'-''^trT,r^im,-D^)W^{C)Tout{m2-D2))A 

(1) 

where Ki = J^^ dTi(mi + ^in,out = 1,75, ••• are meson vertices, and 

Wf^{C) is the Wilson loop with spin insertions, taken along the contour C 
formed by paths {Dz)xy and {Dz)xy, 

W^{C) = PpPAeMw I A^dz^) exp(^ (J^^^F^^dr^ - g H af^F^^dr^). 

ti C "0 

(2) 

The last factor in defines the spin interaction of quark and antiquark. 
The average {Wa)A in (HD can be computed exactly through field correlators 
{F {!)... F{n)) A, and keeping only the lowest one,(F(l)F(2)), which yields 
according to lattice calculation [|15[ accuracy around 1%, |]16| one obtains 

(W^,(C))A^exp(-i[/' ds^,{l) j dsxa{2) + 

+ E /" ^l^dT, r ^''£dT,]{F,^{l)FU2))). (3) 

The Gaussian correlator (F^j,(1)Fao-(2)) = D^i,^\„{l,2) can be rewritten 
identically in terms of two scalar functions D{x) and Di{x) p|, which have 
been computed on the lattice P] to have the exponential form D{x)Di{x) ~ 
exp{—\x\/Tg with the gluon correlation length Tg ?a 0.2 fm. 

Thus the first term in the exponent @) yields the area law {Wo- =) 
exp(— aSmin), with the string tension a 

a = i y D{x)d^x. (4) 

We concentrate now on this confining term in and quote the result for 
the spin-dependent term at the end of this section. 

The important point to be discussed now is the character of dynamics one 
gets from (^, (1). To this end one should compare two characteristic lengths 
(times) Tg and Tq where Tq is the typical period of quark motion, e.g. 
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the classical period of motion along the Coulomb orbit (for heavy quarks) 
or along the orbit in the linear potential (for light quarks). In all cases one 
gets Tg > 1 fm and hence Tg S> Tg. Thus one has a local dynamics however 
relativistic for light quarks, but in any case the excitation of gluonic vacuum 
degrees of freedom can be neglected in the first approximation, so that the 
dynamics is of the potential type. 

As the next step one introduces the einbein variables /ij and z/; the first 
one to transform the proper times Si,Ti into the actual (Euclidean) times 
U = zf. One has 

rt+. roo P 

2/i.(t.) = ^, dsiiD^z^%y = const j Din{ti){D^ z^'^),y. (5) 

The variable v enters in the Gaussian representation of the Nambu-Goto 
form for Smin and its stationary value has the physical meaning of the 
energy density along the string. In case of several strings, as in the baryon 
case or the hybrid case, each piece of string has its own parameter z/*^*-*. 

To get rid of the path integration in (|lD one can go over to the effective 
Hamiltonian using the identity 

G,q{x,y) = {x\eM-HT)\y) (6) 

where T is the evolution parameter corresponding to the hypersurface chosen 



for the Hamiltonian: it is the hyperplane = const in the cm. case |T3 
and Zj^ = const in the light-cone case [^, |T 



The final form of the cm. Hamiltonian (apart from the spin and pertur- 



bative terms to be discussed later) for the qq case is [20 



i=l 



2/i, 2 J 2[/ii(l-C)^ + /i2e + /o'W-C)V 

Here ( = {fii + Jq /3iyd(3) / {fii + yU2 + Jo P^'iP) and Hi and z/(/3) are to be 
found from the stationary point of the Hamiltonian 

-n -n (R) 
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Note that Hq contains as input only mi,m2 and a, where are current 
masses defined at the scale 1 GeV. The further analysis is simplified by the 
observation that for L = one finds u^^^ = ar from (|^) and f^i = \/m^ + p^, 
hence i^o becomes the usual Relativistic Quark Model (RQM) Hamiltonian 

m 



ll^{t = ^)=Y.^ml^^''^ar. (9) 



i=\ 



For large L however one can neglect /ij as compared to v and one gets 

El^2.o^U^\ u^'\(3) = J^ ] (10) 

V vr /i _ 4(^5 _ 1)2 



From (pup one can see that u^^^P) is indeed the energy density along 
the string and jj^^^ is the cm. energy of the quark, which plays the role of 
constituent quark mass as will be seen below. 

To proceed one can use two approximations. First, replace Hq in (|) by 
its eigenvalue Mq, which is accurate within 5% ||22|. Second, approximate 
the L-dependent term in (^, introducing the correction SHgtr namely 

2/2 2 \ 

Ho^Hn + AHstr, ^« = E f + f j + (H) 
and the mass correction due to AHgtr for equal mass case is [|13| 



A,,,,(nL) = (AiJ,,,) = (12) 

where Mq is the eigenvalue of H^] a more accurate approximation is given 
byig 

A„.,„,) . JE1^L±^, ,13, 

But iJo is not the whole story, one should take into account 3 additional 
terms: spin terms in (|^) which produce two types of contributions: self- 
energy correction 



2 l\rr?{i\ An 
Hseif = T.^r^^ A^? = -AK), r^(0) = 1-0.9 (14) 

i=l ^l^i 
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and spin-dependent interaction between quark and antiquark Hspm 
which is entirely described by the field correlators D{x), Di{x), including 
also the one-gluon exchange part present in Di{x). 

Finally one should take into account gluon exchange contributions J^, |11 



which can be divided into the Coulomb part Hcoui = — and Hrad 
including space-like gluon exchanges and perturbative self-energy corrections 
(we shall systematically omit these corrections since they are small for light 
quarks to be discussed below). In addition there are gluon contributions 
which are nondiagonal in number of gluons Ug and quarks (till now only 
the sector = was considered) and therefore mixing meson states with 
hybrids and glueballs ||26|. we call these terms Hmix and refer the reader to 
p6| and the cited there references for more discussion. Assembling all terms 
together one has the following total Hamiltonian in the limit of large and 
small Tg-. 

H = Hq + Hself + Hgpin + Hcoul + H^ad + Hmix- (15) 

We start with Hq = Hr + H string- The eigenvalues Mq of Hr can be given 
with 1% accuracy by |^ 

3 

Ml ^ 8(tL + Anain + -) (16) 

where n is the radial quantum number, n = 0, 1, 2, ... Remarkably Mq ~ 4/io, 
and for L = n = one has /io(0, 0) = 0.35 GeV for a = 0.18 GeV^, and /iq is 
fast increasing with growing n and L. This fact explains that spin interac- 
tions become unimportant beyond L = 0,1,2 since they are proportional to 
dTidT2 ~ 4^i7^'^^i^^2 (see (||) and [^). Thus constituent mass (which is ac- 
tually "constituent energy") fj,Q is "running". The validity of /xq as a socially 
accepted "constituent mass" is confirmed by its numerical value given above, 
the spin splittings of light and heavy-light mesons [29| and by baryon 



magnetic moments expressed directly through fiQ, and being in agreement 
with experimental values . 



The next topic is Regge trajectories in QCD. As it is clear from one 
has the correct asymptotic Regge slope coinciding with the string picture, 
while for small L,L < 4, the approximation (|ll]), (12) holds [23| with almost 



the same slope. The Regge intercept depends strongly on the term Hsei / 
since (0) yields too large value for Mq, e.g. Mo(0,0) ~ 1.4 GeV is almost 
twice the p mass. However the self-energy term (|Hp defined unambigiously 



through a |24] has a negative sign and a magnitude which brings the mass 
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back near the experimental value. Thus one can understand the origin of 
the large negative phenomeno logical constant, which is usually introduced in 
RQM, but it is also rewarding that it is not actually constant, but depends 
on n, L via fio{nL) so that the linear Regge behaviour is preserved. 

To compare with the experiment and disentangle the contribution of spin 
interaction we shall consider the center-of-gravity (cog) masses for each me- 
son multiplet as in Then the masses of all orbital excitations {n = 0) 
can be nicely described by the linear Regge trajectory, which we call the 
Regge L-trajectory with experimental parameters ||2^ 



or 



M^{L) = (1.23 ± 0.02)L + 0.37 ± 0.02(GeV2) (17) 



L = 0.81M2(L) -0.30. (18) 
This is different from the leading p-trajectory 

J = a'jM^iJ) + 0.48, a'j = 0.88GeV^ aj(0) = 0.48 (19) 

since its parameters depend on spin interactions. Now using (^^-jTB]) one 



obtains 23 



L = 0.80M2(L) - 0.34 (20) 

which agrees with experiment (p!8D within 10%, the accuracy being in accor- 
dance with the estimates of the neglected terms in ([T3|), namely H^ad and 
Hmix- In this way one solves the problem of Regge trajectories for orbital ex- 
citations in QCD, thus supporting the foundations o of the wide and fruitful 
activity undertaken by Professor K.A.Ter-Martirosyan and his group to de- 
scribe the high-energy scattering and production processes in the framework 
of the Regge theory. This important contribution was reviewed in . 

We now come to the gluon- containing systems, hybrids and glueballs. 
Referring the reader to the original papers |22[-|Q one can recapitulate the 



main results for the spectrum. In both cases the total Hamiltonian has the 
same form as in ([T3|), however the contribution of corrections differs. 

For glueballs it was argued in [3|] that Hq (pUf ) has the same form, but 



with rrii = and a — >• aadj = ^cr while Hgeif = due to gauge invariance, 
and Hcoui is small due to strong cancellation between tree graps and loop 
corrections [^. Thus glueball masses are expressed through only cr = 0.18 
GeV^ (fixed by meson Regge trajectories) for the center-of-gravity and in 



addition through for the spin splittings. One can see in Table 1 of [34| the 
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theoretical c.o.g. mass values computed in in comparison with lattice 
data. The agreement is striking, especially if one takes into account that in 
theoretical calculations there are no parameters at all -ex was fixed beforehand 
at the same value, as in lattice, and as was neglected altogether for the reason 
stated above. 

We now come to the delicate and very important topic of glueball trajecto- 
ries and especially of the pomeron trajectory. Since the glueball Hamiltonian 
is basically the same as for mesons, one expects that the asymptotic slope of 
all Regge trajectories would be 

a'a{M^ -.00)^ = \a'^. (21) 

One expects that the pomeron trajectory is passing through the states 2+"*" 
(2.29 GeV for a = 0.18 GeV^) and 4++ (around 3.2 GeV) in which case 
the pomeron intercept appears too low ap(0) < 1. A possible way out was 
suggested in [Q, where the intersection of two lowest meson trajectories 
with the pomeron trajectory was introduced yielding the correct intercept 
ap(0) = 1.1. 1.2 for reasonable parameters of trajectory interactions. How- 
ever this topic is far from clarity, and the BFKL perturbative results for ap(0) 
seem to be not stable with respect to inclusion to higher orders and/or non- 
perturbative effects. It seems that the pomeron trajectory depends on both 
perturbative and nonperturbative contributions and their possible interfer- 
ence, and the problem was never considered in that fullness. 

We now coming to the final topic of this talk -hybrids and their role in 
hadron dynamics. We start with the hybrid Hamiltonian and spectrum. This 



topic in the framework of FCM was considered in [^, ^] The Hamiltonian 
Hq for hybrid looks like [H 



i^r^ = !I± + !!± + f^l+J^l+J^ + pI + P^ + a V |r, - r.l + H,,.. (22) 
° 2/ii 2/i2 2 ^ 2/x ^ ' ^ ^ 

Here p^, are Jacobi momenta of the 3-body system, Hs^if is the same 
as for meson, while Hgpin and Hcoui have different structure |^ . 



The main feature of the present approach based on the BPTh, is that 
valence gluon in the hybrid is situated at some point on the string connecting 
quark and antiquark, and the gluon creates a kink on the string so that two 
pieces of the string move independently (however connected at the point 
of gluon). This differs strongly from the flux-tube model where hybrid is 
associated with the string excitation ClS db whole. 
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The difference between two approaches is especially pronounced in the 
case of the hybrid with static quarks where the flux-tube model predicts 
for large R the terms En{R) ~ while in FCM there is another branch 



corresponding to the longitudinal d.o.f. of gluon E{R) ~ p, and 
the lattice results certainly prefer the latter and contradict the flux-tube 
asymptotics. Also at intermediate interquark distances R the spectrum of 
FCM approach is in much better agreement than only another model, see 



[|33| for details and discussion. 

Results for light and heavy exotic 1 ^ hybrids are also given in and 
are in agreement with lattice calculations. Typically an additional gluon in 
the exotic {L = 1) state "weights" 1.2^1.5 GeV for light to heavy quarks, 
while nonexotic gluon {L = 0) brings about 1 GeV to the mass of the total 
qqg system. 

It is understandable now that hybrids play a very special role in QCD, 
namely they describe the excitation of the film - the string world sheet - 
which is between q and q and in particular covers the gluon and quark loops 
appearing in the as renormalization. The "additional gluon" mass Mg ^ 1 
GeV^ cited above enters as a screening mass in the one-loop running as (see 
1^ ||36| for details and derivations 



4% 

"^(^) = — T^m^ (23) 

and this form (in r-space) is in a perfect agreement with a recent calculation 
of as on the lattice i37[| . 

Thus one can say that the perturbative QCD in the IR region is defined 
by the hybrid physics. A similar conclusion can be drawn with respect to 
the DIS diagrams at small to moderate Q^,Q^ ^ 1 GeV^, where familiar 
ladder diagrams with gluon exchanges are replaced by (multi)-hybrid Green's 
functions, and to the hadron-hadron scattering, where e.g. a sequence of 
BFKL ladders is replaced by multihybrid diagrams. The progress in this 
direction is highly desirable. 

All discussion above refers to the leading large limit, which has a 
remarkable accuracy, as it was recently demonstrated on the lattice and 



by comparison of calculated and experimental masses in pj, [M, p3|, B 



We now turn to the 1/Nc effects which as shown in strongly modify 



the masses of radially excited mesons. Typically those mesons have radii 



larger than 1.5 fm and as it was argued in |39] at the interquark distances 



9 



i? > 1.4 fm there appear sea-quark holes in the confining film connecting 



q and q trajectories. It was assumed in that the film with the virtual 
holes establishes the quasistationary state which can be called the "predecay 
state". It is clear that the effective string tension decreases due to holes 
and the confining potential is partly screened at > 1.4 fm. The masses 
of radially excited mesons for L = 0, 1,2,3 and rir = 0, 1,2,3,4 have been 



calculated in with the help of this potential and are shown in Fig.l. 
One can see a remarkable agreement of theoretical masses (black dots) with 
experimental candidates. Another interesting feature is the almost exact 
linearity of trajectories M^(nr, L) vs n^. The effect of the mass decrease due 
to the sea-quark holes is significant (for high it is around 0.5 GeV). 

Concluding the talk I would like to stress the simplicity of the method 
(FCM) which solves at least 4 important problems in the QCD spectrum: 
1) Regge slope 2) Regge intercept 3) the problem of constituent masses 4) 
problem of radial Regge trajectories. 

The author is indebted for useful discussions to A.B.Kaidalov, Yu. S. 
Kalashnikova. During all years of the development of the present method 
the author felt a constant interest, a helpful attention and support from 
K.A.Ter-Martirosyan, which are very gratefully acknowledged. The present 
work was partially supported by RFBR grants 00-02-17836 and 00-15-96786 
and by INTAS grants 00-110 and 00-00366. 
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Figure 1: The spin-averaged meson mass squared M'^{L,n^) vs n = + 1. 
Theoretically calculated masses are depicted by small black dots and are 
connected by straight lines. Other points are experimental candidates. The 
lines from bottom up correspond to L = 0, 1, 2, 3 respectively. 
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